Decoherence and entanglement in a bosonic Josephson junction: Bose-enhanced 

quantum-Zeno control of phase-diffusion 
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We study the effect of decoherence on dynamical phase diffusion in the two-site Bose-Ifubbard 
model. Starting with an odd parity excited coherent state, the initial loss of single particle coherence 
varies from small bound oscillations in the Rabi regime, through hyperbolic depletion in the Joseph- 
son regime, to a Gaussian decay in the Fock regime. The inclusion of local-site noise, measuring the 
relative number difference between the modes, is shown to enhance phase-diffusion. In comparison, 
site-indiscriminate noise measuring the population imbalance between the two quasi-momontum 
modes, slows down the loss of single-particle coherence. Decoherence thus either enhances or sup- 
presses phase-diffusion, depending on the details of system-bath coupling and the overlap of decoher- 
ence pointer states with coUisional-entanglement pointer states. The deceleration of phase-diffusion 
due to the coupling with the environment may be viewed as a many-body quantum-Zeno effect. 
The extended effective decay times in the presence of projective measurement, axe further enhanced 
with increasing number of particles A'", by a bosonic factor of \fN in the Fock regime and N/ log A'' 
in the Josephson regime. 

PACS numbers; 05.30.Jp, 03.65.Xp, 03.75.Mn, 42.50.Xa 



I. INTRODUCTION 

Decoherence and entanglement constitute the last fron- 
tiers of quantum mechanics. The crucial role of decoher- 
ence in the quantum-classical transition, was gradually 
revealed over the past three decades P, 0, |3, Ac- 
cording to this 'environmentally induced superselection' 
scenario 0, 0, the environment interactions with 
an open quantum system, select a preferred basis set of 
so called pointer states that are immune to decoherence. 
Pointer states define realizable observables, thereby re- 
stricting the number of meaningful measurements to a 
handful of classical quantities. For a multipartite system, 
all correlations but those existing in the pointer states, 
arc rapidly destroyed by decoherence, so that out of all 
possible decompositions of the mixed multipartite state 
into a statistical mixture of pure states, only the decom- 
position into pointer states is resilient. Hence, after a 
short time from its onset, decoherence drives the system 
into an incoherent sum of eigenstates of particular clas- 
sical observables. 

The notion of classicality is fundamental to the theory 
of quantum Bose gases Q. In fact, the very definition of 
Bose-Einstein condensation (BEG) is the emergence of 
single-particle coherence poj, which enables the Gross- 
Pitaevskii (GP) classical field description of the interact- 
ing many-particle BEG state. 
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In the above, ?/'(r, t) is the classical condensate order pa- 
rameter, m is particle's mass, Vea;t(r) is an external po- 
tential (e. g. an harmonic trap), and U = iTrh^a/m is 
the two body interaction strength (set by low-energy s- 
wave scattering with characteristic scattering length a). 
Thus, somewhat ironically, the epitomization of quantum 
mechanics on a mesoscopic scale, is in fact a classical, un- 
correlated state. 

At zero temperature, the magnitude of quantum fluc- 
tuations around the condensate state, scales down as 
1/\/N, with increasing particle number N. Since TV in 
dilute BEG experiments is in the range of 10^ — lO'' par- 
ticles, quantum fluctuations are characteristically small 
and the classical field description is an excellent approx- 
imation. However, the initially small value of quantum 
fiuctuations does not necessarily guarantee that they will 
remain small away from equilibrium. In particular, if the 
classical dynamics is unstable, fluctuations grow rapidly 
and become significant on timescales that only grow loga- 
rithmically with N 12, 13, 14, 15]. Single-particle 
coherence is thus dynamically lost due to interparticle 
entanglement caused by the interactions, at a rate which 
depends on the stability of the classical dynamics. For 
large N, the early stages of this process may be viewed as 
the decoherence of the reduced single particle state due 
to its coupling to a bath of quantum correlations [l3l ■ As 
such, collisional dephasing will have its own set of pointer 
states, untouched by the coupling to the 'environment' of 
quantum correlations. 

Goupling the many-body system to an external bath 
(e.g. at finite temperature, where the bath consists of 
the thermal-cloud degrees of freedom) , we are confronted 
with an interesting situation, where the decoherence of 
the system due to its entanglement with the bath, may af- 
fect the loss of single-particle coherence due to the inter- 
nal entanglement of different particles, thereby affecting 
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its classicality. In previous work we found that local site 
noise enhances the collisional loss of single-particle coher- 
ence in a two mode BEC, also known as phase- diffusion 
[11 [il Q E H 113, H H IM Ip. Hence such 
noise drives the system away from classical behavior . 

Here we provide the details of recent work 14] show- 
ing that decoherence can also suppress the loss of single- 
particle coherence in the interacting two-mode system. 
Since inelastic collisions measure the relative number be- 
tween the modes, the pointer states of collisional dephas- 
ing are relative-number (Fock) states. Starting from an 
antisymmetric coherent state in the absence of an ex- 
ternal bath, it dephases in a rate which depends on the 
coupling parameter k = UN/ J, where J is the hopping 
energy between the modes [15|. Subjecting the same 
state to different forms of decoherence, we find that when 
the pointer states of the external bath coincide with the 
relative number states, dephasing is enhanced. On the 
other hand, when the pointer states have a definite num- 
ber between the odd- and even superpositions of the two 
modes, so as to include the initial state, dephasing is 
slowed down. Since the amplification of quantum noise is 
sensitive to the relative phase between the odd- and even 
superpositions of the modes, the suppression of single- 
particle dephasing can be viewed as a quantum Zeno ef- 
fect (QZE) [lllllllllli,!!^,!!!!!!, resulting from the 
continuous projection onto relative number states which 
prevents this phase from taking a definite value. The pre- 
dicted suppression of phase-diffusion is shown to depend 
on the number of particles as N / log N in the 1 < k < N"^ 
Josephson regime and as ^/N in the n > A^-^ Fock regime, 
making the QZE dramatically more pronounced for high 
N . We note that another example of increased coherence 
due to the combined effects of lossy dissipation and in- 
teractions, has been found in the bimodal BEC system 
in the form of a stochastic resonance [13, [Ml . 

The model system is presented in Sec. II and collisional 
dephasing is discussed in Sec. Ill, demonstrating different 
qualitative behavior in three distinct coupling regimes. 
In Sec. IV, we open the system and study the influence of 
two decoherence mechanisms on the phase-diffusion pro- 
cess. Summary and conclusions are presented in Sec. V. 



II. THE TWO-SITE BOSE-HUBBARD MODEL 
AND DYNAMICAL EQUATIONS 

Starting out as a bare-bones toy model [H, [H, [H, [H, 
IbgI . IstI . [sl, m, HO, HH , which retains many of the qual- 
itative features of the many-mode Bose-Hubbard model 
[H [H, Hi, Hil , the two- mode BEC was recently shown 
to be useful in the description of quantum interference 
experiments lip. [ZtI. [isl. l49l| . bosonic Josephson junctions 
[sol Isil Is^ . I53L l54| and quantum tunneling in an array 
of separated double wells [55, 56.,^]. The tight-binding 
condition for such wells is 

\U\n^ = , (2) 



where n is the number density, I = \/h/{muj) is the char- 
acteristic trap size and w is the trap frequency in either 
well. Thus, provided that I 3A^|a|, the system is accu- 
rately described by the quantized Josephson Hamiltonian 

H = (a\a2 + dial^ ^ y ("1 ^ "2) + ^ ("1 ^ "2)^ , 

(3) 

where a.^ and al are the annihilation and creation opera- 
tors respectively, for boson particles in the mode z = 1, 2 
with corresponding particle number operators fii = a]^ai, 
A is a bias potential between the modes, which in the 
following is set to zero, and J is the matrix element cou- 
pling the modes. In the above we have set h = I, and 
eliminated c-number terms proportional to the conserved 
total number of particles N — fii + ri2- 

The quantized Josephson model is easily mapped onto 
a sp in problem by defining the three SU{2) generators 
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which determine the reduced single-particle density ma- 
trix (SPDM) 
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where 1, ctx, cry , <Tz are the identity and Pauli matrices 
respectively, and i,j — 1,2. With these definitions, the 
Hamiltonian ^ assumes the form 



H = -JLx - AL, 
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As shown in Fig [U we account for the effect of two 
types of noise on the two-mode dynamics, by using the 
quantum kinetic Master equation, 
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where p is the A^-particle density operator. The L, noise 
term on the r.h.s. of Eq. © corresponds to the dephas- 
ing of the two local (quasi-coordinate) modes, as may be 
caused by collisions with thermal particles [5§|. The Lx 
noise term is the same for the phase between the odd- and 
even combinations of the site-mode functions (i.e. the 
two quasi- momentum modes). This site- indiscriminate 
noise may be implemented e.g. by a stochastic modu- 
lation of the potential barrier between the sites. Since 
the odd quasi-momentum mode has a node between the 
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Lz noise 

FIG. 1: (a) The (-1,0,0) initial state on the Bloch 
sphere subjected to phase-diffusion, local noise, and site- 
indiscriminate noise. Phase-diffusion and local noise drive 
the Bloch vector towards the Sz-axis, while the Lx noise 
projects it on the Si:-axis. (b) Lz noise distinguishes between 
the left and the right well states, while Lx noise distinguishes 
between the even and the odd quasi-momentum states. 

sites whereas the even quasi-momentum mode does not, 
such perturbation will serve to measure the population 
imbalance between them. Setting A = 0, the equations 
of motion for the angular momentum operators are 



Lx — 



dt 



-U[LyL 



= JL 



dt 



+ LzLyj — V zLx , 

u (l,Lx + - (r, 

</ Ly ^xLz ■ 



^x)Ly , 
(10) 



The classical limit of the two-mode model ([8|) is ob- 
tained by taking the expectation values of Eqs. (|10p and 
using the Gaussian (coherent) approximation 
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which is equivalent to replacing the spin operators Li by 
classical c-numbers. Defining the single-particle Bloch 
vector 
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rescaling time as r = Jt, and substituting into (jlOp . we 
obtain the Bloch equations: 
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where Si = dsi/dr, n = UN/ J is the coupling parameter 
and 7i — Ti/ J are the rescaled decoherence rates. 

For ^x — lz = 0, the mean-field equations depict 
the unitary rotations of the Bloch vector with conserved 
norm |s| — 1. Beyond the mean-field approximation how- 
ever, single-particle coherence will be lost due to entan- 
glement between the particles, even at zero temperature 
where no noise is present. Since we are interested in 
the effect of phase-noise on this entanglement process. 



we must use a quantum formalism and probe its conver- 
gence to the classical (mean-field) limit with and with- 
out noise. A partial account of the dynamics of quantum 
fiuctuations, is provided by the Bogoliubov Backreaction 
(BBR) formalism [H [H |60| , based on the higher-order 
truncation scheme, 

(444) ~ {L,L,){Lk) + {U){LjLk) + {L,U){L,) 

-2(4>(4)(ifc) • (14) 

Using approximation (|14p to truncate the hierarchy of 
dynamical equations for correlation functions, we obtain. 
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-f K [SzAzx + SxAzz) 



-IzAyz - A-lx (Ay^ + SySz) , 

Azx = -Axy - K{SzAyz + SyAzz) - ilz+lx) Azx , 

where the two point normal correlation functions Ay are 
defined as, 
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Finally, a full quantum solution may be obtained nu- 
merically in the Fock representation of mutual , Lz 
eigenstates. 



\Lm) 



1 



y/ {I + my. {I -my. 



l+m 



satisfying. 



L^\l,m) = l{l + l)\l,m) 
Lz \l, m) — m \l, m) , 



|0>, (17) 
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with m — —I, . . . ,0, . . . ,1. The Hamiltonian ([8]), as well 
as the operators producing the decoherence in the master 
equation ([2]), commute with due to particle- number 
conservation. Consequently, the total angular momen- 
tum is fixed at I = N/2. We shall therefore denote the 
Fock states below as |m) = |Z = N/2, m). Using this basis 
set of -I- 1 relative-number states to represent /5, 
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the master equation ^ takes the form, 
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dr 



^Jl{l + 1) - m'(m' + l)p,n.m'+i + VKl + 1) - m'{m' - 
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^Jl{l + 1) - m(m + l)p,„+i,,„' - \//(/ + 1) - m(TO - l)p,„_i 

+ 1) - (m + 1)2)2 _ ^ l)2p„+2.m' 

+ -\/ (/(Z + 1) - (m - 1)2)2 _ _ l)2p„_2,™' 
+ \/ + 1) - (m' + 1)2)2 _ ^ i)2^^_^,^2 
+ + 1) - (m' - 1)2)2 _ _ 1)2^^_^,_ 
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The numerical solution of Eq. (PT|) will be used in the 
following sections to verify analytic predictions based on 
the linearization of the BBR equations ([T5)l . 
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III. PHASE DIFFUSION 

It is clear from the discussion in section II, that the 
classical states of the two-mode system coincide with the 
SU{2) coherent states 
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where < 9 < tt and < cj) < 2tt are the rotation 
angles of the state \m — —I) corresponding to the point 
s = (0, 0, -1) on the SU{2) Bloch sphere. The mean-field 
factorization ([TT]) for these states, is accurate to order 
1/iV, becoming exact for large A^. 

Interactions between particles lead to phase-diffusion 
[11 [il, H, E H [13, H [il H ill 113 . For repulsive 
interactions {U > 0), the collisional loss of single-particle 
coherence due to interparticle entanglement exhibits par- 
ticularly rich dynamics [l^, [H, around the excited 



corresponding to all atoms in the excited quasi- 
momentum mode (and hence to a well-defined relative 
phase (j) = tt). The Bloch vector affiliated with the state 
is So = (—1,0,0). Subjecting it to interactions, one 
obtains different qualitative behavior of phase-diffusion 
in the three distinct interaction regimes [IB]: (a) the 
weak-interaction Rabi regime < k < 1 , (b) the strong- 
interaction Josephson regime 1 < k < N'^, and (c) the 
extreme strong-interaction Fock regime k > N'^ . 

In what follows, we quantitatively characterize the 
early stages of phase-diffusion in the three interaction 
regimes. Separating s{t) = Sq -t- Ss(t) and linearizing the 
mean-field equations (|13p about this stationary classical 
point with "fx — "fz — 0, we obtain 
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resulting in the secular equation A(A2 
the generally complex eigenvalues 
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1 — k) = with 
I (24) 
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and the corresponding nonorthogonal natural modes 

60 = 5sx 

S± — cosQSsy ±isinQSsz, (25) 

where tan O = — k (under this definition, the angle 
Q and its trigonometric functions are real in the Rabi 
regime and imaginary in the Josephson and Fock regimes, 
where they may be replaced by the real angle iQ and 
its real hyperbolic functions). For U > and k < 1 
(i.e. < 8 < 7r/4), all the characteristic frequencies are 
real, indicating the stability of the point Sq. Classical 
fluctuation remain bound as, 

<5+(r) = 5+(0)exp(-zAT), 

(5_(r) = 5_(0)exp(+iAr), (26) 

with A = — k\. Transforming back to Ssi, we obtain 
that fluctuations evolve like the sine and cosine of At, 

Ss^ir) = 6s,i0) (27) 
^^vi''') — '^Si/(0)cos(At) + (5s2(0)tan9sin(AT) , 
Ssz{t) = — (5sy(0) cot sin(AT) + (^52(0) cos(Ar) , 

so that 

(5+(5_ = |(5+|2 = \S^\^ =cos'^edsl + sm^edsl, (28) 



is a constant of motion. In contrast, for 17 > and 
K > 1, imaginary characteristic frequencies appear and 
the point Sq becomes dynamically unstable, 

S+{t) = <5+(0)exp(+AT), 

5-{t) = (5_(0)cxp(-At). (29) 

Consequently we obtain, 

Ss.,{t) = Ss^iO) (30) 
i5sy(r) = (5sy(0) cosh(AT) + (5sz(0) tanh(i6) sinh(AT) , 
Ssz{t) = (5sz(0)cosh(AT) + (5s;,(0)coth(ze)sinh(AT) , 

with tanh(i8) = itanO = ^/ k — 1 and coth(z8) = 
— icot9 = 1/-\/k — 1. Thus, while in the Rabi regime 
the magnitude of (5_|_ and (5_ is conserved separately, only 
the product S^S- is conserved in the Josephson and Fock 
regimes. 

The classical stability analysis of mean-field trajecto- 
ries is invariably associated with the dynamics of quan- 
tum corrections. In order to track down the initial growth 
of these fluctuations, we linearize the BBR equations (fT5|) 
with = Jz — 0, about the point (—1,0,0), to obtain 
the block-diagonal form 
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In addition to the conserved variance A^x , the linearized 
equations (|3ip have the pertinent natural modes. 



small oscillations. 



Ao = cos^ eAyy -f sin-^ eA^^ 
A± = cos^ ©Ayy - sin^ QA^z ± i sin(2e) Aj,^ 



(32) 



with the respective characteristic frequencies 0,2A±. 
Substituting the initial values Ayy = A^z = '2/N, A^y = 
Azx = Ayz = Axx = of quantum fluctuations for the 
coherent state |7r/2,7r), we have 



Ao(r) = 1 , A±(r) = A cos{2e)e^^^\ (33) 



Therefore the variances A± are conserved separately in 
the Rabi regime, whereas the An variances carry out 



Axxir) 

Ayy(T) 
^zzir) 



0, (34) 
2Ao + (A+ + A_) 1 + cos(2e) cos(2Ar) 



4 cos2 e 
2Ao-(A+ + A_) 
4 sin^ 6 



N cos2 e 
1 - cos(2e) cos(2Ar) 

N sin^ e ■ 



Since cos(20) > for < < 7r/4, relative-phase fluc- 
tuations Ayy are reduced at the expense of increased 
relative- number fluctuations A^z- The relative-phase 
squeezing of the excited state \tt/2,7t) is complemen- 
tary to the well-known relative-number squeezing of the 
ground state |7r/2,0) 

In order to relate the dynamics of correlation functions 
in (p4)) to the loss of fringe-visibility in matter-wave in- 
terference experiments, we note that in the absence of 
decoherence (7^ = jz =0), the BBR equations (fT5|) con- 
serve the sum -I- {Axx + Ayy + Azz) /2, which for the 
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FIG. 2: (color online) Absolute value of single particle co- 
herence as a function of rescaled time, starting from the co- 
herent state |7r/2,7r) with (a) k — 0.5, (b) k = 2, and (c) 
K = oo (i.e. UN — 2, J — 0). Solid black, dashed blue, dotted 
red, dash-dotted green, bold solid magenta and bold dashed 
cyan correspond to growing TV in all three figures, with A'' = 
100, 120, 150, 200, 300, 600 in (a), A = 50, 100, 200, 400, 800 in 
(b), and N = 16,36,64, 100, 144, 196 in (c). Gray lines cor- 
respond to the analytic forms of Eq. (|35p in (a), Eq. (|36p in 
(b), and Eq. gl]) in (c). 



initial coherent state is 1 -)- The off-diagonal coher- 
ence of the reduced SPDM is thus, 



|i?i2(r)|^ = 



1 



4 4 
1 cot2(2e) 
4 



2 

sin^(AT), 



(35) 



where L+ = +iLy — d\a2. The fringe- visibility in the 
weak-interaction regime remains high and carries small 
oscillations whose amplitude is inversely proportional to 
the number of particles N. The two modes are phase- 
locked in the Rabi regime. 

Repeating this process for strong interactions k > 1, 
the classical instability is associated with a rapid diver- 
gence from mean-field theory. This deviation is mani- 
fested in the squeezing of the initially coherent state ac- 
cording to Eq. ([33]), by exponential reduction of the A_ 
variance at the expense of increasing A_|_. Both Ayy and 
Azz are amplified hyperbolically and the single-particle 




(b) 



FIG. 3: (color online) Propagation of mean-field equations 
P3p for = 7z =0, with a distribution of initial conditions 
around the unstable point s = (—1,0,0). The initial distri- 
bution is shown in (a). The distribution at r = 2 is shown in 
(b) 



coherence is lost as, 



|i?i2(r)|^ 



coth^(2ie) 



sinh^(Ar) 



(36) 



The phase-diffusion rate is therefore independent of the 
number of particles, but its onset time scales logarith- 
mically with N [l^ [l3| . The growth of purely quantum 
correlations, related to the coUisional entanglement be- 
tween particles in the absence of an external bath, will 
thus vary depending on the stability of the classical dy- 
namics, determined by the value of the coupling param- 
eter K. 

In Fig. [5] we plot the numerically calculated magnitude 
of the single-particle coherence |i?i2|, as a function of the 
rescaled time r, for three values of k, corresponding to the 
different interaction regimes. All the results in this figure 
are obtained for 7^ = 72; = 0. In the Rabi regime (k < 1), 
represented in Fig.[2Ia) by the value k = 0.5, all charac- 
teristic frequencies are real, the relative phase is locked, 
and the fringe-visibility exhibits the anticipated stable 
small oscillations, in excellent agreement with Eq. (|35p . 
The oscillation amplitude decreases reciprocally with in- 
creasing A^ as iV ^ 00, keeping k = UN/ J fixed. There- 
fore in the Rabi regime, small 0(1/A^) initial quantum 
fluctuations remain small compared to the C(l) classical 
mean-field. 

By comparison, correlations grow rapidly in the 
Josephson regime. In Fig. [2jb) we plot the evolution 
of single particle coherence, starting from the state (|22p . 
for K = 2. Since the characteristic frequencies of the lin- 
earized equations become imaginary in this regime, A4. 
fluctuations are amplified at the expense of attenuating 
A_, as described in Eq. ([33)1 . and the initially coher- 
ent state is squeezed. The relation between this squeez- 
ing and the classical dynamical instability is illustrated 



7 




(a) 



(b) 



FIG. 4: Absolute values of the full density matrix entries 
in relative Fock basis, for exact evolution with A'' = 30 and 
K = 2, without noise, (a) Initial coherent state ||- , tt) at r = 0. 
(b) The creation of the cat state at r = 2.3, approximately 
given by (|m = I) + \m — —l))/\/2 (disregarding the relative 
phase). 



in Fig. [3l where the initial quantum propagation is ap- 
proximated by mean-field evolution with a distribution of 
classical initial conditions which replace quantum fluctu- 
ations. While not as accurate as a complete stochastic 
calculation, this crude form of the truncated Wigner ap- 
proach [bI, clearly shows how the classical saddle- 
point leads to the dynamical squeezing of quantum fluc- 
tuations and the amplification of quantum corrections as 
expected from Eq. p6p . Since the deviation from mean- 
field classicality is proportional to sinh^(Ar)/A^, the char- 
acteristic quantum-break time T{,r (the time at which this 
deviation crosses some predetermined threshold of 0(1)), 
only grows logarithmically with N, leading to significant 
phase-diffusion at short times, even for a large particle 
number. We note that this behavior differs from the 
phase-locking of the ground {cf> = 0) coherent state in 
the Josephson regime (l5j . The hyperbolic dependence 
of Eq. p6|) gives an accurate description of the deple- 
tion at short times and a good estimate for the quantum 
break-time. At longer times, depletion becomes suffi- 
ciently large so that the linearized BBR equations (pij) 
are no longer adequate to describe the dynamics. The 
full nonlinear BBR equations however, remain valid and 
trace out the loss of single-particle coherence with good 
accuracy, as shown in Fig. [51 Partial and full revivals are 
observed due to the finite number of {N + 1) phase-space 
dimensions. 

It is also instructive to explore phase-diffusion in the 
Josephson regime beyond the semiclassical squeezing pic- 
ture. In Fig[4]we plot the entries of the iV-particle density 
matrix p with N — 30, for the initial coherent state ([^^ 
and for the evolved state at r — 2.3. Interestingly, in 
the Josephson regime phase-diffusion of the odd |7r/2,7r) 
coherent state with repulsive interactions, results in the 
generation of a macroscopic Schrodinger's cat state (a 
macroscopic superposition of the two Fock states with 
all particles occupying either mode). The same observa- 



tion was made for the even coherent state |7r/2, 0), evolv- 
ing into a macroscopic superposition state under phase 
diffusion with attractive interactions [gI, [H, [g^] ■ 

In the Fock regime, tunneling is effectively turned off 
and the number distribution is locked at its initial Gaus- 
sian form. Time evolution is thus restricted to phase- 
oscillations of the various Fock amplitudes. In this case, 
approximate analytical expressions may be found for the 
loss of multi-shot fringe visibility. Setting J = 0, the time 
evolution of the initial coherent state due to the purely 
collisional Hamiltonian is 



\ip{t)) = exp{-iHt) ^,t:) = cxp{-iULlt) 



TT 

2' 



= exp [—iUrn^t) 



in— — I 



with the binomial coefficients of Eq. (P^ . 

(-1)'+" 



2' 



2' 



21 

I + m 



1/2 



(37) 



(38) 



The expectation value of is thus 



L^) = J2 <„,oC„oe'^(""-"')*(m'|i+|m) (39) 

m' ,7n— — l 
I 



Since relative-number fluctuations for the state |7r/2,7r) 
are Poissonian, we should only consider non- vanishing 
c,„o with m smaller than a few ^/N — \/2l standard 
deviations, so that for large I (i.e. large N) we have 
m,m+l <C I for all relevant m. Therefore, we can replace 

A^/2. Moreover, for large I 



we can replace C(„i_|.i)o 



Cmo, so that. 



le 



iUt 



E 

m— — l 



I 1 2 2iUmt 
\Cm0\ e 



(40) 



Thus, for large N, the time evolution of the single- 
particle coherence is just the Fourier transform of the 
initial number distribution. With the standard large- 
number replacement of the binomial distribution by a 
normal (Gaussian) distribution of the relative number 
difference with width ANr = ^fN, we finally obtain. 



1 



|i?r.(r)| = -exp^-2^^ 



(41) 



with the characteristic dimensionless decay time of 
J/{U\/N) = \/N/k. This behavior is illustrated in 
Fig. IHc), where the gray lines corresponding to Eq. 
(HT|) coincide with the numerical calculations. The same 
dependence of the single-particle coherence decay time 
on (d^/ dNi) i^.^iy f2^Nr (in our case {d^i/dNi)pf.^pf/2 = 
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U and AiVr = vA'') was recently used to demon- 
strate squeezing via interference experiments [2ll |. Sub- 
Poissonian number fluctuations with ATV^. = \A/V / s 
{s being the squeezing parameter) resulted in slower 
relative-coherence loss during a purely coUisional hold pe- 
riod, than would be expected from a coherent state with 
ANr = VN. Gaussian decoherence was also predicted 
for the case of a single spin coupled to an environment of 
N other distinguishable spins [64I . Equation (|^T|) should 
be viewed as a second-quantized version of the same re- 
sult. 

While the Gaussian decay of Eq. (|4T|) is irreversible, 
it is clear that due to the finite dimensionality of the 
system, coherence will revive. From Eq. ([5^ . it is clear 
that Ri2{t + 2-kIU) = i?i2(<) and that i?i2(t + vr/C/) = 
—Ri2{t). Thus, relative-phase coherence with a tt shift 
should be restored after the revival time of — n/U. 
Such collapse and revival was observed in interference 
experiments in optical lattices [i^l , where coherent states 
prepared in the Rabi regime, were allowed to evolve after 
the lattice height was raised, thereby afflicting a sudden 
transition to the Fock regime. 

Summarizing this section, the dynamics of single par- 
ticle coherence starting from the coherent state |7r/2,7r) 
in the absence of noise, depends on the interaction 
strength. In the Rabi regime phase is locked and the 
fringe-visibility does not decay but carries out small os- 
cillations. For stronger interactions the fringe visibility is 
lost on a timescale which increases as the classical limit 
is approached by increasing N while keeping k fixed, as 
logA^ in the Josephson regime and as ^/N in the Fock 
regime. 



IV. DECOHERENCE - QUANTUM ZENO 
SUPPRESSION OF PHASE DIFFUSION 

Having established the different phase-diffusion 
regimes, we proceed to study the effect of noise on this 
process. We repeat the numerical calculations of the 
previous section, using the master equation (|2ip with 
non- vanishing values of the decoherence rates jx and 72 . 
Single-particle coherence is evaluated here using the von 
Neumann entropy of the reduced SPDM, 



S = -TriRlnR) 



1 



■In 



(l + |s|) 



(1- 



(42) 

The effect of local noise induced by collisions with 
thermal-cloud atoms, has been considered in previous 
work [12, [13 ■ From section III, it is evident that the 
pointer states for inelastic collisions within the BEG, are 
also the Lz eigenstates \m). Thus, since both entan- 
glement and decoherence drive the system towards the 
same many-body states (see Fig.[TJa)), such decoherence 
enhances the deviation from classicality. This effect in 
the Rabi regime and in the Josephson regime, is demon- 
strated in Fig. [5] and Fig. [5] respectively, where the dy- 




FIG. 5: (color online) Single-particle von Neumann entropy 
as a function of the rescaled time, for k = 0.5 (weak interac- 
tions). Initial conditions as in Fig. [2] Solid, dash-dotted, 
and dashed lines correspond to numerically exact calcula- 
tions with A'^ = 100, 150, 300 particles respectively. Non- 
Hamiltonian time evolution with 7z = 10~*,7a; — (bold 
blue lines) is compared with the Hamiltonian evolution with 
7z = 7i = (normal red lines) and to mean-field decoherence 
with 7z = 10"^, 7i = (dotted black line). The correspond- 
ing BBR approximations for the Hamiltonian case, are shown 
in gray. 



namics of single-particle coherence in the presence of lo- 
cal noise "fz > 0,"fx = 0, is compared to its evolution 
when 7z = 7a; = (the BBR approximate calculations of 
the decoherence-free, Hamiltonian evolution are shown 
in gray). 

While single-particle coherence is not conserved even 
in the classical limit because 7^ 7^ 0, the deviation from 
mean-field theory in the Rabi regime increases due to 
decoherence (e.g. the minima of the single-particle en- 
tropy oscillations in Fig. [5] lie above the mean-field en- 
tropy). Hence, phase-diffusion and local noise bootstrap 
to give faster dephasing. This trend is even more signif- 
icant in the Josephson regirne (F'ig.JHl). The calculations 
in Fig. [6] differ from Ref . [l2, i only in the choice 
of the initial state, which was taken here to be |7r/2,7r) 
as compared to |0,0) in [13, but since starting from 
the non-stationary |0,0) state drives the system towards 
|7r/2,7r), where most dephasing takes place, the results 
are qualitatively similar. The main conclusion remains 
that for strong interactions (see also bold blue lines in 
Fig. Ilip the quantum breaktime with 7^ > saturates 
to a finite value as compared to its logarithmic growth in 
agreement with Eq. p6p . when jz — 0- 

The saturation of the dephasing time in the Joseph- 
son regime results from the known sensitivity of the 
dynamically-produced macroscopic superposition state 
to local noise [gl, IH, 113] ■ The energy spacing be- 
tween even- and odd cat states is exponentially small 
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FIG. 6: (color online) Single particle von Neumann entropy 
as a function of rescaled time, for k = 2 and 7^: = 0. Ini- 
tial conditions as in Fig. [2] Solid, dashed, dotted, and dash- 
dotted lines correspond to TV = 50, 100, 200, 400 particles re- 
spectively. Non-Hamiltonian time evolution with 72 = 0.1 
(bold blue lines) is compared to the Hamiltonian evolution 
with 7z = (normal red lines) and to the mean-field deco- 
herence with 7^ =0.1 (bold dotted black line). BBR approx- 
imate solutions for 72 = 7a: = 0, are shown in gray. 



0.15 




0.05 



(a) 



(b) 



FIG. 7: Absolute values of the A'^-particle density matrix en- 
tries in relative Fock basis, for exact evolution with A'^ = 30 
and n — 2, with noise, (a) The dephased cat state, created 
by local noise at r = 2.3, for 7z = 0.05, fx = 0, (b) The cat 
state is not formed and the one-body coherence is preserved 
by non-local noise at r = 2.4, for 7^ = 0,"/x = 1. 



in K, so that any small width incurred to these lev- 
els due to the local noise, may couple them and de- 
stroy the macroscopic coherence between the sites. The 
decohered state, shown in Fig. Ela), is close to the 
mixture(|m=/) (m=/| -I- |m= — I) {m— — Z|)/2, correspond- 
ing to a 50% probability of finding the particles in ei- 
ther site [b^, and to the origin s = (0, 0, 0) of the Bloch 
sphere. 

While local noise enhances the dephasing between the 
two modes, a different behavior may be anticipated if 



we consider other forms of decoherence. The multi-shot 
fringe visibility is proportional to the projection of the 
Bloch vector onto the s^Sy plane. Since s lies parallel to 
the Sx axis, we have R12 = Sx/2. From Eq. (j35p and Eq. 
([36)1 we see that for all interaction regimes this projection 
at r < 1/A scales quadratically rather than linearly in 
time, 



|s|=2|i?i2| = 



N 



1 - {u^rf 



(43) 



where u = ^2/iV| coth(2ie)|A. Thus, like the decay 
rate of any quantum-mechanical observable, the phase- 
diffusion rate vanishes at short times. Therefore, the 
rapid phase-diffusion of the initial unstable state in the 
Josephson and Fock reg imes, may be inhibited by a quan- 
tum Zeno effect [2l[2l[2i,[2l,[23[2l[2|]. In its discrete 
form, frequent projective measurements of the relative- 
number between the odd- and even quasi-momentum 
modes, may be carried out by addressing the overlap 
region of the two mode wavefunctions, separating the 
vanishing odd superposition from the non- vanishing even 
combination (see Fig. [T]). If these projections are taken 
at sufficiently short intervals with respect to the Tc = 1/A 
correlation time 5t <^ < they lead to exponential 
decay of the form. 



|s(t)| = [1 - {ujirfY « eH("'5^)'] = pj 



(44) 



with j being the number of successive measurements and 
jSr = r. The effective characteristic break time in the 
presence of decoherence 



1 



Tbr = 



7Vtanh2(2ie) 1 



should be compared in the Josephson regime with the 
decoherence-free breaktime. 



log [7Vtanh^(2ie)/4] 



2A 



(46) 



or in the Fock regime with the r^r = v N / k characteris- 
tic time. Thus, on top of the standard extension of the 
decay time by frequent measurements, manifested by the 
comparison of the denominators in Eq. (|^5)) and Eq. ([^5]) 
with X5t <C 1 , there is a bosonic many-body factor of or- 
der N/ log ^ 1 in the Josephson regime and of order 
^/N in the Fock regime. This factor results from the sta- 
bilization of the initial state which transforms the scaling 
of the breaktime with A^ from logarithmic to linear. The 
many-body QZE will thus be significantly enhanced by 
bosonic stimulation. 

The continuous version of the many-boson QZE 
amounts to the inclusion of the noise term proportional 
to 7a;. Provided that 7^; A > w, we can adiabatically 



(47) 



eliminate A^^ in Eqs. (|3ip and obtain 



|s(t)| = exp ( - — 
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FIG. 8: (color online) Single particle von Neumann entropy as 
a function of rescaled time, for k — 0.5, 7z = 0, and — 100. 
Initial conditions as in Fig. (2] The Hamiltonian dynamics 
with 7a; = (solid red) is compared to the evolution with 
non- vanishing (bold blue lines). Solid, dashed, dotted, and 
dash-dotted bold blue lines correspond to 7^, = 1, 5, 10, 20. 
Gray circles, squares, triangles, and diamonds respectively, 
portray the quantum-Zeno prediction of Eq. (|47p . substituted 
into Eq. (|42p . Analytical values coincide with numerical re- 
sults for 7 ^ A. 



xlO 




X 



FIG. 9: (color online) Single particle von Neumann entropy 
as a function of rescaled time, for « = 0.5 and 7z = 0. Initial 
conditions as in Fig. (2] Solid, dash-dotted, and dashed lines 
correspond to A'' = 100, 150, 300 particles respectively. Non- 
Hamiltonian time evolution with = 10 (bold blue lines) is 
compared to the Hamiltonian evolution with 72, = (normal 
red lines). Anticipated QZE behavior for the same particle 
numbers, is depicted by gray circles, squares, and triangles 
respectively. 



so that the characteristic depletion time is 

T6,. = % > - > Ur- (48) 

The QZE suppressed phase-difFusion is illustrated in 
the numerical results of Fig. [S] and Fig. [S] for the Rabi 
regime and in Fig. (TU] and Fig. [TT] for the Josephson 
regime. In Fig. [8] and Fig. [TOl the initial evolution of 
the single particle entropy at large jx, is shown to ap- 
proach the predicted value, obtained by substitution of 
Eq. (|T7|) into Eq. Both the weak- interaction os- 

cillations and the strong-interaction phase-difFusion are 
replaced by an exponential (i.e. initially linear) decay, 
with a rate proportional to l/^N"/^), because they are 
piecewise constructed from similar short-time building 
blocks. The continuous projection onto the axis, may 
thus either destabilize the stable excited state in the Rabi 
regime (Fig. [5]), or halt its unstable decay in the Joseph- 
son regime (Fig. fTO)) . From a quantum-to-classical transi- 
tion perspective, when applied to the strongly-interacting 
system, this form of decoherence, which is tantamount 
to the continuous measurement of the relative-number 
in the rotated basis set of quasi-momentum states, will 
prolong the characteristic depletion time, thus restoring 
mean-field classicality. 

Bose enhancement of the QZE is illustrated in Fig. [3] 
and Fig. llOi showing for weak- and strong coupling re- 
spectively, the QZE with varying particle numbers. Since 
the decay rates and the initial slopes are indeed de- 
pendent on 1/A'^, it is evident that the transition from 
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FIG. 10: (color online) Single particle von Neumann entropy 
as a function of rescaled time, for k = 2, 72 = 0, and A'^ — 100. 
Initial conditions as in Fig. (2] The Hamiltonian dynamics 
with 7;c = (solid red) is compared to the evolution with 
non- vanishing "fx (bold blue lines). Solid, dashed, dotted, 
and dash-dotted lines correspond to 7^, — 1,5,10,20. Gray 
circles, squares, triangles, and diamonds respectively, portray 
the quantum-Zeno behavior of Eq. (|47|l for the same values of 
7a; , substituted into Eq. (|42|l . Similar to the weak-interactions 
case (Fig. |8} good agreement with the QZE prediction is at- 
tained for 7 ^ A. A clear transition from initially quadratic 
to initially linear depletion is observed. 
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FIG. 11: (color online) Single particle von Neumann entropy 
as a function of rescaled time, for k = 2 and 7^ — 0. Ini- 
tial conditions as in Fig. [2] Solid, dashed, dotted, and dash- 
dotted lines correspond to = 50, 100, 200, 400 particles re- 
spectively. Non-Hamiltonian time evolution with = 10 
(bold blue lines) is compared to the Hamiltonian evolution 
with "fx = (normal red lines). Anticipated QZE behavior 
is depicted by gray circles, squares, triangles, and diamonds. 
QZE is more pronounced for higher particle number. 

logarithmically growing breaktimes followed by an A^- 
independent dephasing rate (normal red lines in Fig. Ill|) 
to 1/N dependent decay rates (bold blue lines), implies 
the enhancement of the QZE by the expected factor, with 
stronger suppression of phase-diffusion for larger values 
of A^. 

Looking at the full density matrix (Fig. Efb)), we see 
that the site-indiscriminate noise protects the initial co- 
herent state, retaining a localized distribution somewhat 
smeared along the main and secondary diagonals. Thus, 
in contrast to the local noise which figuratively speak- 
ing 'kills the cat', the projection onto aborts the 
cat's very birth. Also by contrast, while even weak site- 
specific noise suffices to destroy the fragile coherence of 
the macroscopic superposition, the parity/momentum- 
specific noise needs to be sufficiently strong to induce 
the QZE effect. 

To summarize our results, the quantum breaktime 
growth with increasing N (keeping k — 2 fixed in the 
Josephson regime), is plotted in Fig. [121 In the absence of 
decoherence, we obtain the expected logA^ convergence 
depicted by the bold black line (BBR) and circles (ex- 
act numerics, in good agreement with BBR). When local 
on-site noise (72 > 0, bold red lines and symbols) is ap- 
plied, this convergence is slowed down as shown earlier 
in Fig. [6] and the quantum breaktime saturates to a finite 
value which decreases logarithmically with increasing 7^ 
[T^ [Tsj . In comparison, when the noise is focused in 
the overlap region between the single-particle site-mode 
functions so as to measure the number-difference between 
their nonlocal odd- and even superposition (7^; > 0, nor- 
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FIG. 12: (color online) Time at which S reaches 0.2, starting 
from the coherent state |7r/2,7r), as a function of the particle 
number A^, for n = 2. The bold solid black line, calculated 
from the BBR equations (|15|l . depicts the log TV dependence 
for 72, = 7z = 0, in good agreement with the exact numerical 
results (black circles). Bold red lines correspond to BBR cal- 
culations with 7a; = and 7z = 5 x 10~* (dashed), 5 x 10~^ 
(dotted), 5 X 10~^ (dash-dotted with inverted triangles denot- 
ing exact numerical results) , all showing slower than logarith- 
mic convergence with saturation to a finite breaktime at large 
A^. The QZE transition to linear convergence is illustrated by 
the normal blue lines (BBR) and symbols (exact numerics), 
corresponding to 7z = and 7a; =0.1 (dashed, diamonds), 
0.4 (dotted, squares), and 1 (dash-dotted, triangles). 



mal blue lines and symbols), we observe faster than loga- 
rithmic convergence, approaching the QZE linear depen- 
dence on N for large 7^;. 



V. CONCLUSIONS 

Using the two-site Bose- Hubbard model, we have stud- 
ied the collisional phase-diffusion of the excited quasi- 
momentum number state with all particles in the odd 
superposition of the two modes. The loss of single- 
particle coherence is associated with classical stability, 
with bound oscillations in the stable Rabi regime and 
rapid depletion for strong interactions, in agreement with 
[13, fT5 |. The strong-coupling instability is presumably 
associated with the rapid heating observed in the merg- 
ing of two condensates with a tt relative-phase, on an 
atom chip [4§]. Linearizing the BBR equations, we ob- 
tained exact forms which describe the initial BEC de- 
pletion with good accuracy. The loss of coherence, at- 
tributed to the bose-amplification of initial spontaneous 
emission noise from the excited odd-superposition state 
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to the even ground state, was shown to be quadratic (i.e. 
non-exponential) in t within the — 1/A correlation 
time. 

In order to study the effect of decohercnce on phase- 
diffusion, we have considered two types of noise, cor- 
responding to the measurement of the relative-particle 
number in two different basis sets. In the double-well 
realization, relative population imbalance between the 
site-modes can be measured by local noise on either site, 
whereas the relative number difference between the quasi- 
momentum modes (i.e. the even- and odd superpositions 
of the site-modes) can be determined by probing the over- 
lap region in between the sites. The first type of deco- 
herence was shown to accelerate depletion, in accordance 
with previous results [H, By contrast, odd-even rel- 
ative number measurement, results in the suppression of 
coUisional dephasing. This somewhat counterintuitive 
behavior, where one form of dephasing effectively sup- 
presses another, is essentially a quantum Zeno effect, 
wherein the decay of an unstable state is suppressed by 
its continuous observation. However, unlike the classic 
QZE suppression of spontaneous decay [IJ, the bosonic 
enhancement of the initial quantum noise implies that 
in subjecting the system to decoherence, a transition is 
made between iV-independent decay rates, to rates that 
scale linearly with N. Consequently, the deceleration 
of phase-diffusion becomes more effective as the number 



of particles is increased by a N/ log N factor ( Josephson 
regime) or by a \/N factor (Fock regime), which may 
become very significant for realistic particle numbers in 
dilute gas BEC experiments. 

In addition to the realization of the two-mode dy- 
namics in the split ting and mergiiig of atomic BECs 

Ei 113, m, H M,M, m, H, m, M, 113. similar 

behavior should be obtained in a wide range of many- 
boson quantum systems, which exhibit classical dynam- 
ical instabilities. One such example is the stimulated 
dissociation of a molecular BEC into bosonic fragments 
[H, IM IZO, [Zil, Izi [Zi, [Zi. where the short-time dy- 
namics of the atom-molecule number difference depict 
a similar amplification of spontaneously-emitted atom- 
pairs, and therefore adhere to the same non-exponential 
behavior [H, [t^ . Frequent or continuous observation of 
the relative atom-molecule number, will hence suppress 
this stimulated process, with similar bosonic enhance- 
ment factors. 
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